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Abstract 

The statistical behavior of a nonlinear system described by a mapping with phase rotation is 
studied. We use the Kolmogorov-Chapman equations for the multi-time probability distribution 
functions for investigation of dynamics under the external noise perturbations. We find a stationary 
solution in the long-time limit as a power series around a state with complete phase randomiza- 
tion ("phase mixing"). The Ornstein-Uhlenbeck and Kubo- Andersen models of noise statistics are 
considered; the conditions of convergence of the power expansions are established. 
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1 Introduction 

Noise is a key governing factor in behavior of dissipativc nonlinear systems. Strictly speaking, any 
purely deterministic description appears to be inconsistent, because dissipative processes are connected 
to random kinetic motion (thermal processes, spontaneous emission, radioactive decay, etc.), so that 
they are stochastic in nature. Nevertheless, deterministic models provide an fair description of averaged 
motion for systems with stable trajectories. In the opposite case of unstable motion, random fluctua- 
tions may be amplified, so that validity of non-stochastic descriptions becomes questionable. The latter 
case is of particular importance, because many effects, like noise-induced transitions ]lj and determin- 
istic diffusion j^, are inherent in stochastic nonlinear systems. It is also worth to mention current 
attempts to reconsider the background of statistical mechanics and to demonstrate that deterministic 
chaos underlies the origin of macroscopic irreversibility (|]-||. 

In this work we consider an evolution of a nonlinear noisy discrete dynamical system. We concentrate 
on evolutionary mappings with a controllable mixing, exemplified by the Ikeda mapping 0. By this is 
meant that the maximal Lyapunov exponent L is an increasing function of a control parameter A and the 
case L,A> 1 corresponds to chaotic motion with intense phase mixing The main idea of our work 
is to use A^ 1 as a small parameter and to analyze power series around a "basic" solution corresponding 
to the limiting case A — > oo. It is important that in the limit A — » oo the original nonlinear problem 
reduces to a linear one ||. Analytical description of such limit transition is important by itself and it 
may be useful for validity check of computer simulations. As it was shown in |], [l^] for the case of 
intense phase mixing, the stationary solutions are associated with certain multifractal measures. 
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We present some rigorous results for a nonlinear circuit with delayed feedback, associated with 
the Ikcda mapping, which served as a model for radiation field dynamics in an optical ring cavity 
with a nonlinear medium H, O]. It can be viewed as a particular case of a wide class of mappings 
with controllable mixing. For instance, reduction of ordinary differential equations to mappings (for 
systems under the action of piecewise constant or delta-shaped forces) often leads to the mappings with 
" nonlinear phase rotations" . In this situation the intense mixing occurs if a small displacement of a 
phase point leads to a non small "rotation" represented via an element of the orthogonal or unitary 
group. The 3D mapping for a nuclear magnetization system fl2| , [l3| , the 2D mapping for parametrically 
driven spin waves in ferromagnetic [fl4f , the 2D mapping for a periodically kicked dissipative oscillator 
[[iJJ belong to this class of mappings. 

The random dynamics of the system under consideration is determined by external noise statistics. 
Due to the delayed feedback, the noise contributions from the different times form a superposition. For 
this reason it is important to preassign time correlation properties of the external noise unambiguously. 
Here we consider two widely used models of external noise: the Ornstein-Uhlenbeck jlj and Kubo- 
Andersen |l6| random processes. Both of them may be thought of as an outcome of the Brownian 
motion in the thermostat with continuous or discrete assembly of states. 

In the Section 2 we introduce a nonlinear model, write the generalized Kolmogorov-Chapman [ fl7| 
equations for the multi-time distribution functions and derive the corresponding equations for their 
Fourier transforms. In the next Section, we examine the case of intense phase mixing A — > 00 leading 
to reduced equations; the stationary solution is found by applying the iterative method. In the Section 
4, we find stationary solutions for the original equations in the form of power series. In the Section 5, 
we prove a convergence of series solution and obtain the convergence conditions for this series. In the 
Section 6, we discuss a decisive role of the "rapid" Gaussian noise. 



2 A basic model 

Consider a circuit with a nonlinear element (NE) and a delayed feedback. We assume that the transfor- 
mation of a slowly varying complex valued amplitude X consists of the phase shift (j> — > + A |A| 2 + 0o, 
and the dissipative contraction \X\ — » K \X\, k < 1. The signal travels from the NE through the feedback 
loop (with a round-trip time T d ) and reaches the mixing device, where it adds to the external signal 
1 + £ (t). Assuming £ (t) being a random process with zero mean value (£ (t)) = 0, we can write the 
stochastic difference equation of motion 

X(t) = £(t)+F(X(t-T d )) (2.1) 
= £ (t) + 1 + kX (t — T d ) exp {zA \X (t - T d )\ 2 + z6 } . 

may be rewritten in the form of two-dimensional (complex) mappings family 

X N+1 = £ N + F(X N ), (2.2) 

where X N = X (t + NT d ) and £jy = £ (t + (N + 1) T d ), < t < T d . We suppose that the external 
noise consist of two components - one is the Ornstein-Uhlenbeck random process with small correlation 
time: r cor <C T d (hereafter referred to as the "rapid" Gaussian noise); the Fourier transform of its 

distribution function Pq (rj) is given by A (U) — exp ( — Q \U\ 2 /4 ) . A second noise component is assumed 



Equation (2.1 
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to be the Ornstein-Uhlenbeck or Kubo- Andersen Markovian random process with arbitrary r cor and the 
transition density functions: 

(i) for the Ornstein-Uhlenbeck process 

W fori, [r]) = {nR (1 - C 2 )}" 1 exp (- |c - vCr\ 2 /R (i - C 2 )) , (2.3) 

(ii) for the Kubo- Andersen process 

W (£, v, M) - CrS {2) (£ - V) + q (0 (1 - Cr) , (2.4) 

K 

where ( r = exp (— r/r cor ) , = (£ — £fc) and (a;) = 5 (Res) 5 (Inix) is the 2Z? (complex) 

fc=i 

delta- function. The presence of the separate "rapid" Gaussian noise component in the external signal 
is a decisive element in the convergence proof in Section ^. 

Note that the Kubo- Andersen random process (a special case of the kangaroo process |Q) is a 
stepwise constant process describing random jumps between complex values k = 1,2, ...,Ko (with 
appearance probabilities pk); switching times are uniformly and independently distributed along the 
time axis. 

Under the above-mentioned assumptions we can write the generalized Kolmogorov-Chapman equa- 
tion & ft% for the multi-time distribution functions in the form 



P N+1 ((X s ) , C) = J dYd(,Q ((X s ) , (Y.) , C, Pn ((Y s ) , , (2.5) 

where 

/n 
dYdfjd^---dC n Y[P G (vk) (2.6) 

fc=0 

x<f< 2 > (X k - & - % - F (F fe )) W (&+!,&, [r fc+1 - r fc ]) . 



In Eq.(2.5) we used a short-hand notation 



P N ((X s ) ,£) = P (X [NT d + r ] , X x [NT d + n] , X„ [AT d + r n ] , £ [iVT rf ]) , (2.7) 

where JQ, i = 0, 1, ...,n, are the signal amplitudes at NTd + U, respectively, and £ is the amplitude of 
the second noise component at NTd- We also require = To < t± < ... < t„ + i = Td and use a notation 
dA = dA dAi...dA n . 

The corresponding equation for the distribution function Fourier transform 



*jv ((U s ) ,n) = J dXdt P N {{X s ),i) exp j-*Re (j2 XkU k + } . 



(2.8) 



is given by 

9n+i((U.),0)= /.A- ( /tfTiir,i.irj.!)(i! vj; A ((Vj.^i. (2.!)) 

where 
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T ((U s ) , (V s ) ,Q,6)=Z ((U s ) ,Q,e)f[A (U k ) a (U k ,V k ) , 

fe=0 

a(U,V) = (2ttY 2 J dYexp{iRe(YV* -F (Y )U*)}, 

/n 
dn x ...dn n J[ h (fifc+i.flfc - u k , [ Tk+1 - r k }) , 

k=0 

h (n, r, [r]) = (2^)- 2 /" dtdn w (£, 77, M) ex P {zRc (7?r* - eo*)} . 



(2.10) 
(2.11) 
(2.12) 
(2.13) 

Subscripts replacement N, N + 1 — > si reduces Eqs ■ ( |2 . 5| ) and ( |2.9| ) to the equations for the stationary 
distributions: 



P st ((X s ) , = J dYd^Q ((X s ) , (Y s ) , C, Pst ((Y s ) , . 

((ig , fi) = [dvder ((u s ) , (v s ) , n, 0) * st ((\g , 0) 



3 Intense phase mixing case 

The expression ( |2.1l| ) may be presented in the form || 

a = a + Act, 

where 

ct (U, V) = (2tt |F|) _1 e - lRec/ (5( 2 ) (|V| - « |t/|) , Act ({/, V) = e -* Rec/ /3 ([/, 1/) , 



(2.14) 
(2.15) 



(3.16) 
(3.17) 



1 i 

IS — 1 



i /kIuvi 



x cos < ^ ( 6 + arg 



\V\' + K 2 \U\ 2 (l/+l)7T 

4Az; + 2 



(3.18) 



In the limit A — > oo one can neglect the last term in Eq.( 3.16| ), in other words, replace ct by ct in ( [2.10 ) 
(corresponding stationary solutions of Eqs. (2. 14), (2.15) will be labeled by superscript ^°- ) ). Consider the 
special case of the Ornstein-Uhlenbeck random process in detail; the Kubo- Andersen process is analizcd 
the same way. Using the explicit form of the function ( [2.1S ) for the Ornstein-Uhlenbeck process: 



H (n, T, [r]) = exp |-1 |tt| 2 (1 - C r 2 ) | <5 (2) (T - QCr) . 



(3.19) 



one can reduce Eq.(2.15) to 



Noise Transformation in Nonlinear System 



5 



(3.20) 



where 



and 



/n 
dVd9 f ((V s ) , 9) Z((U S ) , U n+1 ,9) 11 A (U k ) a (U k , V k ) 
i. — n 



71 

n [ A (Vk)> 



iReu k 



k=0 



n+l 



$ ((U s ) , U n+1 ) ( / ( (sC/.e*-) , 

J=0 



We) 



(3.21) 



$((U s ),U n+1 )=exp{-j 



n+l 



k=l 



n+l r 



u, 



C 2 



c 2 



(3.22) 



— p -T k /T cor 



(-)(,) /(-) 



The angular brackets in Eq.( |3.2l| ) denote an averaging over the phase angles <pk, k = 0, 1, ...,n. Define 
a linear operator L acting on the vector [(U s ) , fi] = [Uq, U n , fi] as follows 



L(0o,...,<£»)[(tf«),fi] = 



( K ^ e ^),^Cr fc C/ fc +Cr„ +1 « 



fe=0 



(3.23) 



It is easy to show that 



i ■ ■ ■ i Vnp 



<f> np )}[(U s ),Sl} = 



K m U s exp J i J2 K\ ) > £ + 



fc=0 



(3.24) 



where 



q=0 I p=l J 

and 0fcp are arbitrary angles: < 4>k p < 27r, I < p < m, < k < n (the operators in the product in 
Eq.( [3.24 ) are enumerated from right to left). One can show that 



m — 1 



\d km \ < U £ K'C^- 1- * < Cn ' m (max{ K ,C T „ +1 }) 

9=0 

As consequence, we have 



m— 1 m — >oo 



0. 



(3.25) 
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[(f/ s ),O] = [(Q),0] 



(3.26) 



and 



lim K m /((C/ s ),r!) = /((0),0)- lim K m l 



(3.27) 



for a continuous function / (in the r.h.s. of Eq.(3.27 ) the operator acts to the right onto the real 
constant equal to unity). With help of the expression (3.27), using Eq.(3.2C) in the iteration procedure 



and taking into account the normalization condition 



(3.28) 



*i?((0),0) 



dXd£P£ ) ((X s ),Z) = l, (3.29) 
we can write a solution of Eq.( 3.2(j| ) in a simple form 

¥$P((E7») > n)= lim K m l. (3.30) 

m — >oo 

Note that the solution ( 3.30| ) is the asymptotically stable " fixed point" of the evolutionary mapping 



(o) 

N+l 



(o) 

jV ' 



(3.31) 



so that it is independent of initial conditions. This reflects decay of correlations and "memory loss" in 
the system. The transformation Pjy — * Pn+i defined by Eq.(pTq) preserves the norm for arbitrary A. 



Consequently, the condition |/| < 1 implies 



K/ 



< 



1, so that absolute values of J3.27 ) and (3.30) arc 



bounded from above. The operator K performs transformations including dilation \U S \ — » k|J7 s | and 
averaging over the n— dimensional torus. It is remarkable that we will have Eq. (3.31) in the place of 
Eq.([2.9]) if we substitute the linear stochastic difference equation 



X (t) = £ (f) + 1 + kX (t - T d ) exp {i~6 {t)\ 



(3.32) 



instead of the nonlinear Eq^^); here 9 (t) is the random process with in depen dent and uniformly 
distributed values. As it was shown in our work p0| , the solutions of a type (3.30) are associated with 
certain integrals over multifractals. 

Using statistical moments approach, one may arrive at the same results in other way Define the 
stationary moments of distributions as follows: 



k=0 



where multiple derivatives 



k=0 



Qr k Qh 

dU* k rk dU l k " 



f((u s ),n) 



(3.33) 



(3.34) 



(u s )=o,n=o 
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are computed at the origin, and 7 = (p,q, (r s ,l s )) denotes the multi-index. It follows from Eq. fl3.21| ) 

that -D 7 (K/ J is a linear combination of D > (/) for 7 < I7I (here we denote I7I = p+q + X) ( r k + Ik))- 

^ ' k=0 
Introduce a vector m (/) of values D 1 (/), |7| < |7 mo:r |, and choose special ordering of its components: 
provided that either I72I < I71I, or I72I = |-yi I with p 2 + q2 > Pi + qi, D~ n (/) precedes D l2 (/): 



m (/) = (..., D 7l (/) , —,D~f2 (/)>-)• 

It is easy to see that 

m(K/) =am(/) + b, b = b • / ((0) , 0) , 
where a is a triangular matrix. The diagonal elements 



r k lk K 



rk+h 



(3.35) 
(3.36) 

(3.37) 



fe=0 



with S p9 = e~ Td / Tcor + (l — e~ Td / Tcor ) S p oS q o for the Kubo-Andersen process and S p? = e ~ T<i ( p+9 )/ rc ° r 
for the Ornstein-Uhlenbeck process, are real and < a 77 < 1, so that both a and 1 — a are nonsingular 
matrices. Furthermore, a = TJT -1 , where detT ^ and J is the Jordan form of the matrix a (each 
diagonal element of J also is positive and less than unity). In this case a m = TJ m T~ 1 and 

lim fm(K"7) -(l-a)- 1 b -/((0),0)} 

m — too L V / J 

= lim a™ f m (/) - (1 - a)" 1 b ■ / ((0) , 0)} = 0. (3.38) 
Matrix relations ( 3.36| ) for various 7 ma£C present another equivalent form of the condition (3.27). 



4 Series expansion of exact solution 



In the general case of arbitrary values of A we can rewrite Eq. (|2.15| ) as 



(4.39) 



where e is an auxiliary parameter. We are seeking a solution of Eq.( f4.39| ) in the form of formal power 
series in e mM: 



U) 

St 



with normalization conditions 



3=0 

^ t ((0),0) = *i? ) ((0),0) = l 
^H ) ((0),0) = Q, j>0. 



(4.40) 

(4.41) 
(4.42) 

Substituting the expansion (4.40) into Eq.( |4.39| ) and equating coefficients of like powers of e, we obtain 
Eq.flDol) and 
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Repeated application of the last relation yields 



m— 1 
p=0 



(4.43) 



(4.44) 



where m is an arbitrary positive integer. In the limit m — > oo in Eq.(4.44), using Eqs.(3.27), (4.42) and 
its corollary: 

(4.45) 



lim K m ¥£ =0, j > 0, 



one can find a recursive relation 



p=0 



(4.46) 



Finally, combination of recursions (4.46) with different j and substitution into Eq.(4.40) yields 

oo ( / oo \ ^ ^ 

,(°) xVjJ I \^ k-p I c I ,t,(o). (4.47) 



In the case A 3> 1 the term with in Eq.(4.47) is O (A 2j ); note that one must set e = 1 in the final 
expression. 

5 A proof of convergence 



In this section we find a condition of existence of the expansion ( [f.47| ). Consider the associated equation 
for the "fixed point" \I> st : 



for the operator: 



= A# st , 

a/ = *2 ) + |Ek4s/. 



^p=o 



The operator K is defined by Eq.(3.21) and 

Sf ((U s ) , SI) = J] A (U k ) / dVdO f ((V s ) , 0) Z (([/,) , 0, 6) a ((U s ) , (V.)) 

/ n 



where 



k=0 



fc=0 v ; v ; (r s ) 



i=0 i=fc+l 



(5.48) 
(5.49) 

(5.50) 
(5.51) 
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The inner summation in (5.51) is taken over all the permutations 



1 2 

ro n r 2 



n 

r„. 



For k = n the right-hand product of factors a in (5.51) is equal to unity. We are going to prove that 
A is the contraction operator in an appropriate complete metric space. The condition of contraction is 
given by 



A/ 2 - A/i 



< a con ||/ 2 - /i|| , 



(5.52) 



According to the contraction mapping principle, Eq.(5.48) has the unique solution if the condition ( |5.52| ) 
is satisfied. From Eq.( 5.49| ) follows 



Af 2 - Ah = \J2K" S(/ 2 -/i) 



(5.53) 



.p=o 



and the operator in the right hand side of (5.53) is linear. In this case the condition ( 5.52| ) takes place 
simultaneously with the limiting condition for the operator norm 



def 



sup 

f*0 



£k4s/ 



< Oicon < 1. 



(5.54) 



The ex pans ion ( 4.47 ) is convergent under the condition ( |5.54 ). On the other hand, ( 4.47 ) is the solution 
of Eq.(4.39) resulting from the iterative procedure. 
In this work we use the norm 



ll/l 



def 



sup \\f((u s ),n)\Y[^(u k ) 



(5.55) 



fc=0 



with 



mI') = \MU)}- 1/2 =exp{±Q\U\ 2 



(properties of this norm are discussed in the Appendix). To prove the inequalities (5. 52), (5. 54) and 
make estimate of a co „, we have to make a number of auxiliary estimations. Using Eq.(3.19) and the 
similar expression for the Kubo- Andersen process: 



h (n, r, [r]) = Cr<5 (2) (n - r) + (1 - Cr) x (n) -5 (2) (r) 



(5.56) 



where 



X(n) =^p fc exp{-iRefefi*)} 
fc=i 

(see Eq.( ^.4[ )), one can find for both noise statistics 

d0H(Sl,9-U,[r]) f(U,6) 



sup 

nec 



< sup 1/(17,11)1, 
nec 



(5.57) 
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sup 

nec 



dO Z((U s ),n,6) f((U s ),6) 



< sup|/((E/ s ),n)|; 
nec 



(5.58) 



/ (U, 51) in ( 5.57 ) and / ((U s ) , f2) in ( 5.58 ) are arbitrary bounded complex functions. Next, using ( 3.21 ), 
we can obtain 



sup 
nec 



n 3-1 ; 

K j /((^),0) < TT TT A(^ m (7 fe )(sup \f((^U s e^),n) 
1 fe =om=o \^ ec 



(5.59) 



Besides, using (3. IS) and the inequalities \ J V (z)\ < \z/2\ 1 ' for real values of z and integer positive v 
[19|; k u ^ +1 > < k? +1 , v~ v ~ x < 1 for z/ > 1, j > 0, K < 1, we can find the upper bounds 



2ttA 



|t/v| 

~4A~ 



< 



2ttA 



exp 



|t/| 2 + l^| 2 ' 
8A 



(5.60) 



Finally, the inequalities ( 5.59p , ( ^5 . 60 ) are combined to produce the contraction condition (5.52) with 



Of 



1 — K 



4k 



QA-1 



71+1 



(5.61) 



The parameter a con is less than unity for any chosen n if QX is sufficiently large. Under this condition, 
the solution (4.47) of Eq. (2.15) exists as the unique "fixed point" of the contraction operator (5.49). 
It is possible to prove convergence of (4.47) in another way, using the technique of a majorizing series, 
but it's very cumbersome and we omit it here. 



6 Conclusions and discussion 



We have found a stationary solution to the generalized Kolmogorov-Chapman equation for the multi- 
time distribution functions in the form of a convergent series expansion (4.47). Our approach is based 
on successive application of two iterative procedures. The final expressions were obtained for the case of 
the nonlinear Ikeda mapping and certain fluctuation statistics. This result may be used for calculation of 
the approximate expressions for probability distributions, statistical moments and correlation functions. 

As simple example consider a case with the "rapid" Gaussian noise component only. In this case 
Eq.( [2.15' ) for the Fourier transform of the one-time distribution function takes the form: 



4< st (U) = A (U) / dV V st (V) a (U, V) 



(6.62) 



Using the general expression ( 4.47 ), we can write approximate first-order solution as follows: 



(to = (£/)+£ k4 s* <°> (to = 



iReu 



A(E0{»»oo (k\U\) 



+ / vdvA (v) Voc M ft {U, v) + ^ A (k p U) Vp ^ (k \U\)fo (k p \U\ ,«)]}, 
{ \ p=i I 



(6.63) 
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where 



01 (U,v) = 2iz(0(U,ve i<t ') exp {-rocos0}) (0) , (6.64) 
2 (u,v) = 2tt(0 (ue # \«e^ 3 )exp{-i(ucos0i + vcos0 2 )}) (0i , , (6.65) 



A— 1 



m (u) = TT A (k 3 u) Jo (k 3 u) , T) = 1, 77^= lim (6.66) 

In general case the resulting formulae are cumbersome. Nevertheless, they provide a basis for develop- 
ment of computational procedures (in particular, using symbolic processors). The detailed study of this 
applications is beyond the scope of this article. 

Consider a specific role of the "rapid" Gaussian noise in our method. In order to simplify consider- 
ation, we suppose that this is the only noise component present in the system. To analize influence of 
noise, we replace a (U, V) by A (U) a (U, V) in Eq.(|j|). Since A (U) = exp (-\Q |?7| 4 ) , the Fourier 



spectrum becomes truncated in the region \U\ 3> 1. It ensures the convergence of the expansion (4.47) 
On the other hand, as it was shown in this replacement leads to a "coarse-grained" dynamics, 
because fluctuations coarsen the fine-scale structure of a chaotic attractor. If Q\ is sufficiently large, 



the condition of contraction (5.52) may be satisfied for any n, so that corresponding expansions (4.40), 
(4.47) converge. However, we can consider the limits A — > 00 and Q — > separately. The limits do 
not commute, so that we should set Q — > after using A — > 00. It seems plausible that the correct 
description of the random behavior in the limit A — > 00 may be reached only under the action of an 
infinitesimal noise. 

It is instructive to draw a parallel between the above manipulations and the well-known formalism 
of quantum mechanics. In order to obtain the retarded (causal) solutions of the Schrodinger equation, 
one have to write the equation containing the additional source term: 

i> € + i (h + V) V, = -e (A - *) , (6.67) 

where $ is a solution for the " nonperturbed" system = 0^ , and further to set e — * after the 

volume grows infinitely (this is the boundary condition of Gell-Mann and Goldberger j2(| in the form 
proposed in [2lj [; here h — 0). The same approach appears to be natural to formulate one well-known 
version of the so-called "master equation" formalism in statistical physics. Following the technique of 
pl[ , the final equation of motion for the quasi-equilibrium statistical operator may be derived from the 
Liouwille equation (completed with the source term) by substituting e — ► after the thermodynamic 
limit is reached. Thus, in order to obtain a physically significant irreversible solution, one must add the 
infinitesimal dissipative term (in other words, to break the symmetry with respect to t — > —t). 

Both above-mentioned procedures lead to a " coarsened" description, presented by " coarsened" equa- 
tions. It is worth to note that dissipative processes accompanied by fluctuations bring into existence a 
strange attractor the same time destroying a fine structure of the attractor at A — > 00 (for an example 
of the attractor in our model, see p3| ) . 



Appendix 

In this work we use a nonstandard metrics, so it is necessary to check the completeness of the respective 
metric space. Define a norm and a distance function on a set of functions of n complex arguments 
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x = (xi,X2, ■ ■■,x n ) <E C™ as follows: 



||/|| 00 = SUp |/(x)|, Poo(/l,/ 2 ) = ||/l-/ 2 || t 



xeC" 



(A.68) 



Using the standard approach (see, for example, |2j|, Theorem 1.23), we can prove 

(i) if fk (x), fc G N, is a Cauchy sequence, the pointwise and uniform convergence takes place: 



(Vx G C") (Ve > 0) (3N £ ) (Vfc > N e ) : \f k (x) - / (x)| < e; 

(ii) if (x) is continuous for each k G N, / (x) is continuous as well. 
In addition, we can prove: 

(iii) if all fk (x) are uniformly bounded 



(A.69) 



(3M) (Vfc G N) (Vx G C" ) : |/ fe (x)| < M, 

then |/ (x)| < M. 

It is convenient to prove (iii) by reduction to contradiction. Assume that 

(3x GC"):|/(x )|=M + e, e > 0. 



(A.70) 



(A.71) 



Choosing e < e and k s ufficie ntly large to have (xo) — / (xq)\ < e, we arrive at a contradiction to 
the original assumption ( A.71 ) : 



|/ £ (xo)| >||/(x )M/fc(xo)-/(xo)|| >M + e-s>M. 



(A.72) 



It means that the continuous uniformly bounded functions, equipped with the norm and metric (A.68), 
form the complete metric space. 

It follows that the space C* C C (C" +2 ) of continuous functions, satisfying 



|/((t/ s ),0)| <M\{^ l {U k ), 



(A.73) 



fe=0 



with the norm ( 5.55| ) and the distance pc* (/ii J2) = ||/i — /2II) is the complete metric space (the 
condition ( A.73 ) is rewritten as ||/|| < M). 

Finally, note that / £ C* implies Af G C* . Choosing a class of functions with ||/|| < M, where 
M is sufficiently large: M > (1 — a C on) _1 (the functional space C* may be expanded by replacement 
M — > M , M > M, if required), we may find that 



\\Af\\ < 



J2* p )sf 

\p=o J 



< 1 + a con 11/11 < 1 + a con M < M. 



(A.74) 
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